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Resummation of Cactus Diagrams in Lattice QCD, to all Orders 
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We show how to perform a resummation, to all orders in perturbation theory, of a certain class of gauge invariant 
tadpole-like diagrams in Lattice QCD. These diagrams are often largely responsible for lattice artifacts. Our 
resummation leads to an improved perturbative expansion. Applied to a number of cases of interest, e.g. the lattice 
renormalization of some two-fermion operators, this expansion yields results remarkably close to corresponding 
nonperturbative estimates. We consider in our study both the Wilson and the clover action for fermions. 



1. INTRODUCTION 

Ever since the earliest days of lattice field the- 
ory, one problem present in most numerical sim- 
ulations has been the calculation of corrections 
induced by renormalization on Monte Carlo re- 
sults. Several methods have been used to address 
it. Perturbation theory provides in principle a 
methodical means of calculating renormalization 
functions, operator mixing coefficients, etc. Its 
drawbacks lie in its asymptotic nature, and in 
that it is a formidable task on the lattice, which 
places severe limitations on the order to which 
it can be carried out; indeed, at present, exact 
calculations in perturbative lattice QCD reach 
only two loops (for 2-point diagrams) and 
three loops (for vacuum diagrams) [Q. In re- 
cent years there have been considerable efforts 
in refining the perturbative computations intro- 
ducing recipes motivated by mean field and tad- 
pole resummation arguments J^] . Various nonper- 
turbative, numerical approaches to renormaliza- 
tion functions have also been devised and there 
has been recent progress both in their range of 
applicability and in their precision Non- 
perturbative methods, such as those of Refs. 
H are in general preferable to approximations 
based on perturbative calculations, due to their 
better controlled systematic errors {O(a) against 
O(<?0 ))■ However, improved perturbative esti- 
mates are still quite useful. They indeed provide 
important consistency checks. Further, in those 



cases where nonperturbative methods are difficult 
to implement, perturbative methods remain the 
only source of quantitative information. 




Figure 1: A cactus 

We present an improvement of lattice perturba- 
tion theory, which results from a resummation to 
all orders of a certain class of diagrams, dubbed 
"cactus" diagrams |) 10 1. Briefly stated, these are 
tadpole diagrams which become disconnected if 
any one of their vertices is removed (see Figure 1). 
Our original motivation was the well known ob- 
servation of "tadpole dominance" in lattice per- 
turbation theory. This observation must clearly 
be taken with a grain of salt: One-sided inclu- 
sion of tadpoles can ruin desirable partial can- 
cellations between tadpole and non-tadpole dia- 
grams; worse, their contribution is gauge depen- 
dent. The class of terms we propose to resum 
circumvents the latter objection since, as we shall 
see, it is gauge invariant; it also overcomes the for- 
mer objection in known cases. Given its diagram- 
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matic nature, our method leads to results which 
can be systematically improved, order by order. 
Applied to a number of cases of interest, the cac- 
tus resummation yields remarkable improvements 
when compared with the available nonperturba- 
tive estimates. As regards numerical comparison 
with other improvement schemes, such as boosted 
perturbation theory || , our method fares equally 
well on all the cases studied. 

2. THE CALCULATION 

Consider the standard Wilson action for 
SU(N) lattice gauge fields: 

5 = iE Retr ( 1 - c/ ^) 

U°^ v is the usual product of link variables around 
a plaquette in the fi— v plane with origin at x ; in 
standard notation: 



U 



x^u 



J-goA^,^ igoA^ +IJ ,^ -ig a A mJrUtli -ig Q A x 



where A x>fl = A%^T a . By the Baker- Campbell- 
Hausdorff (BCH) formula we have: 

Kp» = ex P { igo Fi% + 0(g 2 )} 

— A 4- A i — A , —A 

The diagrams which we propose to resum to all 
orders will be cactus diagrams made of vertices 
containing Fxjiii ■ Let us see how such diagrams 
will dress the gluon propagator; we write: 

= + — •— + -*-•- + ..■ 

where the one-particle irreducible piece is given 
by the recursive equation: 




The presence of only Fxjiv m the vertices ensures 
that the longitudinal parts of all propagators can- 
cel, so that the effect of dressing is the same in 
all covariant gauges. Thus, we can write: 



• = w(g ) ■ 

Indeed, the dressed propagator will become a 
multiple of the bare transverse one, where the 
factor w(go) will depend on go and N, but not 
on the momentum. The above diagrammatic re- 
lations can be turned into an algebraic equation 
for w(go) that can be easily solved numerically. 

Similar to the dressing of the propagator, one 
dresses the action vertices, operator insertions, 
etc. Moreover these calculations can be ex- 



tended 1 10 1 to the case of the clover improved 
action formulation of lattice QCD |11|, which is 
widely used in numerical simulations in order to 
reduce scaling corrections. 

Cactus resummation may be applied cither to 
bare quantities or to quantities which have been 
calculated to a given order in perturbation the- 
ory; thus contributions which are not included in 
the resummation can be reintroduced in a sys- 
tematic manner. 

3. APPLICATIONS OF CACTUS RE- 
SUMMATION 

The resummation of the cactus diagrams can 
be applied to the calculation of the renormaliza- 
tion of some lattice operators. Approximate ex- 
pressions of the lattice renormalizations are ob- 
tained by dressing the corresponding one-loop 
calculations. 

As first example |^] we consider the renormal- 
ization Z(go) of the topological charge density lat- 
tice operator 



Q(x) 



2 4 x 32tt 2 



±4 



t r { XJ x . v Ux , p a } 



using the Wilson action and for SU(3). Z(go) 
is a finite function of go, going to one in the 
limit go — > 0, and that is much smaller than 
one in the region go — 1, where Monte Carlo 
simulations using the Wilson action are actually 
performed. In perturbation theory one has || 
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Z(go) = 1 - 0.908#j + 0{g$), leading to the es- 
timate Z(go = 1) ~ 0.092. The large one-loop 
coefBcient suggests already that perturbation the- 
ory can hardly provide an acceptable estimate of 
Z(go) for go ~ 1 without some kind of resum- 
mation. Indeed the heating-method |l3| , which 
does not rely on perturbation theory, gives the 
estimate Z(g = 1) = 0.19(1). Cactus dressing of 
the one-loop expression leads to 

Z { gl) » [l-w(g )} x [l-w(go) - 0.491«? 2 ] 

Since w(g = 1) ~ 0.250, one finds Z(go = 1) — 
0.193, in agreement with the heating-method re- 
sult. A similar estimate is obtained by "mean- 
field" improving § the one-loop calculation: Z ~ 
P 2 (1 - 0Mlg%) = 0.209, where P ~ 0.5937 is 
the plaquette and <?p = gft/P a new coupling. 
Further confirmation || of the validity of the cac- 
tus resummation has been achieved for SU(2). 

One can also apply cactus resummation to 
the lattice renormalization of fermionic opera- 
tors. In particular we consider the renormaliza- 
tions Zv,Za of isovector fermionic currents us- 
ing the clover action jl(J . Their one-loop result is 
known Q, Z V , A = 1 + g 2 zv,A(csw) + ■ ■ ■ Cactus 
dressing amounts to: 




i_i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ~ 

0.2 0.4 0.6 0.8 1 

Si 



In the above figures we show results for Zy 
and Za (from Ref. ||), coming from numeri- 
cal simulations (filled circles, fitted by a solid 
line), bare perturbation theory (dotted lines) and 
"mean field improved" perturbation theory || 
(crosses) . The superimposed dashed lines are our 
results from cactus dressing. As this comparison 
shows, the dressed one-loop expressions consti- 
tute a remarkable improvement with respect to 
the simple one-loop calculation. 
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